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Abstract 

Muon spin relaxation (/xSR) and nuclear magnetic resonance (NMR) experiments revealed that 
the spin singlet state with an excitation gap of ~200 K is realized from S = 1/2 ND4 tetrahedral 
clusters in a cluster Mott insulator GaNt^Ss- The intercluster cooperative phenomenon to the 
singlet state at Ts = 32 K is triggered by intracluster Jahn- Teller type structural instability de- 
veloped from ~3Ts. Referring to the lattice symmetry, the formation of Nbs octamer (Nb4-Nb4 
bond) is suggested. 

PACS numbers: 75.30.Kz, 75.50.-y, 76.75.+i, 76.60.-k 
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In any magnets, if spin is well denned at the atomic site, the classical spin degree of 
freedom must be compensated on approaching K, which is usually achieved by magnetic 
order. Another route via electron-lattice coupling is quenching of spin itself by formation 
of a spin singlet state. It is of interest to extend this concept to cluster compounds with 
mixed-valent magnetic clusters (more than one atom shares unpaired electrons within the 
cluster unit) embedded in a matrix crystal. When the cluster unit has half-integer spin, 
particularly, since intracluster interactions never compensate all the spin, we expect long- 
range magnetic order due to minor intercluster interactions. In this article, we show another 
possibility, formation of a new bound state between clusters, just like dimerization in the 
one-dimensional (ID) S = 1/2 spin chain with the electron-lattice coupling jj. 

As a category of mixed-valent magnetic cluster compounds, ternary calcogenides AB 4 X 8 
(A = Ga, Al, Ge; B = V, Mo, Nb, Ta; X = S, Se) with the cubic GaMo 4 Ss structure 
(space group F43m) {2], 3, 4] are of particular interest. In the compounds, cubic (B 4 X 4 ) n+ 
and tetrahedral (AX 4 ) n ~ ions are weakly coupled in a NaCl manner, resulting in hopping 
conduction between the clusters. GaNb 4 Ss with (Nb 4 S 4 ) 5+ and (GaS 4 ) 5_ ions has attracted 
attention as a sort of Mott insulator, where 4d electrons are localized at Nb 4 clusters due to 
electron correlations, and also as one of pressure-induced superconductors |5(. In molecular 
orbital schemes, the tetrahedral Nb 4 unit in (Nb 4 S 4 ) 5+ shares seven 4d electrons, resulting 
in one unpaired electron at the highest occupied molecular orbital (HOMO) as shown in Fig. 
[T] jsj. In other words, the Nb 4 tetramer has a local moment with S = 1/2. Accordingly, the 
magnetic susceptibility x obeys the Curie- Weiss law with the effective moment fi c g = 1.73 
/ie/f-u. [fj. The negative Weiss constant 6 = —298 K suggests antiferromagnetic (AF) 
coupling among clusters. \ shows an abrupt drop at Tg = 32 K, which is accompanied by a 
structural transition to a tetragonal state (P42im) [6| possibly due to Jahn- Teller (JT) type 
instability; the triply degenerate HOMO splits to singlet and doublet states (see Fig. [[]). 
Concerning the magnetic ground state below Tg, Pocha et al. [5| commented no additional 
reflections in powder neutron diffraction profiles, ruled out long-range AF order, and then 
ascribed the absence of long-range order to topological spin frustration. On the other hand, 
in a more recent paper, Jacob et al. |6( proposed non-collinear AF order on the basis of their 
crystallographic analysis and band structure calculations. Being inconsistent with S = 1/2 
inherent to the cluster, our muon spin relaxation (/iSR) and nuclear magnetic resonance 
(NMR) experiments clearly indicate presence of a spin excitation gap below Tg, namely 
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compensation of spin to the S = state. Referring to the lattice symmetry [5J, we discuss 
possible Nbs octamer formation from two ND4 tetramers. 

All polycrystalline samples were prepared by solid state reaction from pure elements 
sealed in evacuated quartz tubes. Final heat treatment was performed at 900° C for 1 day 
after several heat treatments and intermediate regrinds. \i measured to check sample quality, 
was the same as literature data [3j, |5|, |6| (our result is shown in Fig. 0) Zero-field (ZF) and 
longitudinal-field (LF) /iSR measurements were made at the RIKEN-RAL Muon Facility at 
the Rutherford- Appleton Laboratory in the UK using a pulsed positive surface muon beam 
at 4.5-100 K under LF fields of 0-0.4 T. NMR experiments for 69 - 71 Ga (nuclear spin 1 = 3/2 
for both) and 93 Nb (/ = 9/2) were performed with a standard phase coherent spectrometer. 
Most NMR spectra were obtained by Fourier transform of the spin-echo or free induction 
decay signal in a fixed magnetic field of 7.0071 T. The nuclear spin-lattice relaxation rate 
1/Ti for 69 ' 71 Ga was obtained by inversion and saturation recovery methods. 93 Nb-Ti is too 
short to be estimated reliably. 

To determine the magnetic ground state, we first show results of /iSR. Since NMR is 
affected by both magnetic and electric hyperfine interactions, the analysis of powder-pattern 
spectra is often non-straightforward. In contrast, /xSR probes only the magnetic field, and 
has an advantage in monitoring the appearance of internal fields. Typical examples of ZF- 
/iSR spectra are shown in the inset of Fig. [2j Even at the lowest temperature (T) , no evidence 
of magnetic order such as muon spin precession was observed. Above Ts, the relaxation is 
dominated by Gaussian-like depolarization. The commonly used damped Kubo-Toyabe 
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FIG. 1: Proposed molecular orbital schemes of the Nb4 unit above and below Ts [5(. To interpret 
isotropic and anisotropic hyperfine fields at the Nb site, schematic spin density distributions are 
shown for HOMO. 
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(KT) function P^t) = exp(-Xt)Gf T (A,t) with Gf T {A,t) = § + f(l - AH 2 ) exp(-±A 2 t 2 ) 
was fit to the data, where A/7 M is the width of the static field distribution (7^: muon 
gyromagnetic ratio), and A is the damping rate associated with an additional relaxation 
process mostly due to electron spin fluctuations. Above Tg, by treating A and A as free fitting 
parameters, we obtained a nearly T-independent value of A ~ 0.1 yus -1 , which corresponds 
to a tiny field of A/7 M ~ 1 G, and is ascribed to dipolar fields coming from randomly 
oriented nuclear spins. With decreasing T passing through T s , due to an appearance of fast 
relaxation, the above function tunes to be inappropriate. Taking account of the symmetry 
lowering in the crystal, i.e., possible separation of the muon stopping site, we applied P^{t) = 
(C — 1) exp(— At)G^ T (A, t) + Cexp(— Afast^) to fit low-T data, where C is a constant. In 
addition, associated with appearance of the fast relaxation, the damping part (dynamic 
fields) tends to mask the Gaussian part (small static fields), resulting in unreliable estimation 
of A below ~30 K. Then we assumed T independence of A to analyze low-T data. Thus 
estimated A and Af ast are plotted in Fig. [2] {C gradually increases with lowering T and reaches 
~0.2 at the lowest T). A is enhanced strongly below Ts and approaches a constant value 
at the lowest T. This behavior is commonly seen for the system which condenses into a 
spin singlet state, and now understood as sporadic dynamics due to spin excitations in a 
singlet sea [3, Q]- In LF-/xSR experiments below T s , strong muon-spin depolarization was 
observed in the whole range of LF up to 0.4 T. This observation indicates the presence of 
dynamic components of electronic spins and supports the absence of large static field due 
to static magnetic order. Hence we conclude that the ground state is neither long-range nor 
short-range magnetic order. 

Next let us see NMR data. Examples of 93 Nb-NMR lines measured above Tg are shown 
in the inset of Fig. [31 where only the center line (the m = —1/2 <-» 1/2 transition) is shown. 
Nearby the sharp center line, unresolved satellites due to quadrupole interaction, spreading 
out in several hundred kHz, were observed; the Nb atom occupies an axially symmetric site 
(16e) in the high-T cubic state. With decreasing T, the resonance frequency once decreases 
until ~80 K and increases again. The T dependence of the magnetic hyperfine shift 93 K, 
estimated from the peak position, is shown in Fig. [31 A high-T part of 93 K scales well with x 
llOj. I 93 ^! exhibits a characteristic broad hump at ~80 K, deviating from a Curie- Weiss-like 

n 

T dependence at high T [9|. Furthermore it should be noted that the lineshape turns to be 
strongly asymmetric in spite of the reduction of | 93 if|. The origin of the characteristic hump 
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FIG. 2: (color online). The T dependence of the ZF muon spin relaxation rate A. The inset 
shows typical examples of ZF-/iSR spectra (at 4.5, 25 and 100 K). Solid curves indicate the fit 
by the damped KT function. Below T$, the asymmetry is fit together with another fast damping 
component. 

and anisotropy of 93 K will be discussed later. Below Ts, the 93 Nb spectrum turns to be broad 
and complex. This is mainly due to separation of the Nb site to three crystallographically 
inequivalent sites (8/ and two 4e) and appearance of large anisotropy in 93 K. Here we do 
not discuss further the low-T 93 Nb spectrum, because the analysis is not straightforward 
but complicated. 

Figure |U(a) shows 71 Ga-NMR lines above Ts. A single symmetric line accords with the 
cubic symmetry of the Ga site (4a) above Ts. With decreasing T, the line is broadened. We 
confirmed that the linewidth is magnetic in origin; the ratio of 69 Ga and 71 Ga linewidths 
is T-independent and close to m ^/ 71 ^ = 0.7870 (7: gyromagnetic ratio) but different from 
69 Q/ 71 Q = 1-60 (Q: quadrupole moment). The 71 Ga hyperfme shift, K, estimated from 
the peak position, is also plotted in Fig. [3] [llj. 71 K is in reasonable proportion to x at least 



5 




FIG. 3: (color online). T dependences of 93 Nb and 71 Ga hyperfine shifts, 93 K (solid circles) and 



llfl. The 



71 K (open circles). 93 K is shown only above Tg. 71 K below Tg corresponds to 71 ifj 
solid curve indicates the bulk susceptibility The broken curve represents exp(— Eg/k^T) with 
Eg/k-s — 220 K. The inset shows 93 Nb-NMR lines (only the center line) above Tg, measured under 
an external field of 7.0071 T. The broken straight line indicates the zero shift position. 



above Tg [l0(] in contrast to the case of 93 K. Below Tg, 71 K drops suddenly, and approaches 
a constant value at low T, proving that the low-T upturn in \ is extrinsic. In addition, 
another component extended over 2 MHz appears on the foot of the center line as shown 
in Fig. 11(b). This bottom component is most reasonably assigned to quadrupole satellites. 
According to the lattice symmetry [6|, a single Ga site (nonaxial 2c) exists only even below 
Tg. The unresolved satellites are reasonable if considerable asymmetry appears in 71 K; note 
that the center line looks asymmetric. Focusing on the bottom part of the spectrum, the T 
variation is shown in Fig. H](c). Note that the quadrupole broadening already appears well 
above Tg in spite of the nominal cubic symmetry at the Ga site in the high T state. The 
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FIG. 4: (color online). 71 Ga spectra measured above Ts (a) and below Ts (b), and magnified 
bottom structures above Ts (c). The data were collected in a field of 7.0071 T and plotted with 
respect to peak positions. The spectrum of (b) was obtained by summing Fourier transforms of 
the spin-echo signal accumulated at different frequencies. 

result indicates that local environments around the Ga site already gets distorted well above 
T a . 

Figure [5]( a) shows T dependences of 1/Ti for 69 ' 71 Ga. The recovery of nuclear magne- 
tization fits well with a single exponential function well above Ts, with two exponential 
functions expected for the I = 3/2 nucleus with quadrupole interaction just above Ts, and 
with neither the single exponential nor the 1 = 3/2 function below Ts. These facts are, 
as a whole, consistent with the observation in the static spectrum. The last fact indicates 
the distribution of 1/Ti, which is attributable to the appearance of hyperfme anisotropy. 
Below Ts, 1/Ti was estimated by applying a stretched exponential function exp(— (t/Ti) 13 ) 



(obtained f3 is nearly T independent of ~0.6). With decreasing T, 71 (1/Ti) once decreases to 
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FIG. 5: (color online), (a) T dependences of 1/T\ measured for 69 Ga (open circles) and Ga 
(closed circles). Reliable 69 (1/Ti) was not obtained below T§ due to overlap of 93 Nb and 69 Ga 
signals. The solid curve represents exp(— Eg/k^T) with Eg/k-Q ~ 190 K (b) T dependences of 
quadrupole (open squares) and magnetic (solid squares) components of 71 (1/Ti). 

reach a minimum at ~50 K, increases again to take a sharp peak at Ts and drops abruptly 
below Ts- 

Measurements of 1/Ti for different isotopes 69 ' 71 Ga give valuable information on spin and 
lattice dynamics. As seen in Fig. Ufa), magnitudes of 69 (1/Ti) and 71 (1/Ti) are reversed 
at ~60 K. Because of m 1 / n 1 < 1 and 69 Q/ 71 Q > 1, this implies that both magnetic and 
electric quadrupole fluctuations contribute to 1/Ti, and furthermore depend differently on 
T. 1/Ti for the nucleus with / and Q is expressed as the superposition of magnetic and 
quadrupole terms, 1/Ti = (1/Ti) 7 + (1/Ti)q. Figure Mb) shows T dependences of separated 
contributions, estimated by solving simultaneous equations for 69 (1/T!) and 71 (1/Ti) with 
use of relations (1/Ti) 7 oc 7 2 and {1/T\)q oc Q 2 . The magnetic component decreases 
monotonically, indicating no enhancement towards Ts. For the system with dynamically 
fluctuating local moments, we expect T- independent 1/Ti. The gradual decrease of 1/Ti 
is possibly ascribed to the development of spin-singlet correlations. On the other hand, 
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interestingly, the quadrupole contribution shows a critical behavior on approaching T§. This 
enhancement indicates that structural fluctuations develop from ~100 K, of the order of 3T$. 

ZF-/iSR experiments detected no additional static field below T$. NMR spectral analyses 
are consistent with this observation. These results prove unambiguously that the ground 
state is not static magnetic order. This is in good contrast to the AF state in GeV 4 S§ 
established due to intercluster interactions [151 ] . As seen in Fig. [3j 71 K drops rapidly below 



T s , which is well fit by exp{-E g /k B T) with E g /k B ~ 220 K. As in Fig. [5(a), 71 {1/T 1 ) also 
shows a similar T dependence with E g /k B ~ 190 K. Nearly the same estimation of the gap 
~200 K from both K and 1/T-y indicates that the gap presents in spin excitations. Hence 
we conclude that the ground state of GaNb 4 S 8 is a S = spin-singlet state. The T and field 
dependences of /iSR-A also support this conclusion. 

Characteristic precursory phenomena of the transition are summarized as (i) at the Nb 
site, the isotropic term of the local susceptibility xioc looks suppressed from ~3T S , and 
instead, large anisotropy appears, and (ii) around the Ga site, in contrast to moderate 
magnetic fluctuations, both static and dynamic structural fluctuations are enhanced from 
~3Ts. As shown in Fig. [TJ spin density distributions of HOMO are expected to be highly 
symmetric and asymmetric in the high-T and low-T states, respectively (similar discussions 
have been made on a 3d analog GaV 4 Ss [3, Q)- On this basis, the above results are 
interpreted as follows; intracluster reconfiguration of Nb-4d cluster orbitals (i.e., mixing 
of the low-T electronic state) already starts from ~3Ts, enhances intercluster structural 
instability, and finally induces cooperative long-range symmetry lowering at Ts, which is 
coupled with the spin degree of freedom. 

As the mechanism of the spin-singlet formation, spin compensation via intercluster inter- 
actions should be considered, because each Nb 4 cluster unit has one unpaired electron with 
S = 1/2. The crystallographically unique Nb 4 cluster even below Ts js] rules out charge 



disproportionation among clusters. According to the low-T lattice symmetry [6j, all Nb 4 
tetrahedra are elongated along one of the threefold axes (111) of the cubic lattice. As a 
result, two Nb 4 tetrahedra come closer slightly along the tetragonal [110] or [110] direction. 
Hence, if electron-lattice coupling is appreciable, it is reasonable to expect binding of two 
Nb 4 tetramers to a Nbg octamer (Nb 4 -Nb 4 bond), just like the dimerization of the ID spin 
chain. We speculate that the transition at Ts is related with the characteristic electronic 
state of this material as a cluster Mott insulator [5, E3, 1^], where the electron transfer is 
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suppressed due to electron correlation, but is easily recovered by perturbation. Simultane- 
ously, as mentioned above, we should note that this transition is driven by the local JT-type 
instability within the cluster. In this sense, the transition may be classified as the orbitally 
driven Peierls state [16J, although intercluster spin-singlet correlation may also be developed 
from high-T of the order of spin gap. As a whole, the spin singlet state of GaNJ^Sg is classi- 
fied as one of exotic states found in frustrated lattices near metal-insulator boundaries; the 
crystal structure of GaNh^Ss is closely related to the highly frustrated pyrochlore lattice, 
and simultaneously Nb 4 units arrange as a frustrated fee lattice. 

In conclusion, /iSR and NMR experiments revealed that the ground state of GaM^Ss 
is a nonmagnetic spin singlet with an excitation gap of ~200 K. Taking into account the 
lattice symmetry below Tg, formation of a Nb§ octamer from two Nb 4 tetramers is suggested 
for the first time. The transition is driven by the intracluster JT-type instability developed 
from high T of the order of 3Ts. 
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"Novel States of Matter Induced by Frustration", and Grant-in-Aid for Young Scientists 
(Start-up) from MEXT, Japan. 
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